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I. INTRODUCTION 
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It is well known experimental fact that identical pions, 
which are produced in heavy ion collisions, being bosons 
show Bose-Einstein correlations. These correlations result 
from the quantum mechanical interference in the correspond- 
ing symmetrized n-particle wave function (where n denotes 
the number of produced identical pions). They contain large 
amount of information about the statistical properties of the 
momentum and configuration space distribution of the system, 
and thus provide a potentially very useful method to probe the 
geometry of the hadronizing source (see, for example, and 
references therein). 

According to common understanding we are not able to 
determine which pion is emitted from which position in the 
source, so we are required by Bose statistics to add amplitudes 
for all possible alternate histories. In general, symmetrized 
wave function for n^-pion state can be written in the follow- 
ing way Jl]: 
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where a(j) denotes the / element of a permutation of the 
sequence {1,2, ...,n K } and the sum over a denotes therefore 
the sum over all n K l permutations in this sequence (depen- 
dence on positions of points of detection, which will vanish 
when calculating probabilities, was neglected). Here r's de- 
note the points of production of secondaries. Because in the 
experiment one observes only momenta of produced secon- 
daries these r's must be somehow get rid of. This is so far 
always done by integrating over {r,} with some assumed dis- 
tribution p({rj}), which is assumed to be factorizable and ex- 
pressed by product of independent single particle distributions 
P({ r j}) = YljP( r j) IB Eh- As result one gets the following 
probability of the n^-pion state, 
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expressed by permanent of the matrix 1 14>,' ; ' 1 1, where 

$«y = / e ife - p >p(W)rf 4 {r}. (3) 

For usually considered n% = 2 case one recovers well known 
classical expression for the probability of detecting two pions 



Unfortunately the execution time of direct computation of the 
permanent, eq.(Q~|i, grows exponentially with n n one has to de- 
vise some special methods like Metropolis procedure inves- 
tigated in J2j,j3|] or von Neumann accepting-rejection method 
proposed in Jj] . 

The Metropolis procedure proposed in J2|] uses the stan- 
dard Monte Carlo technique due to Metropolis. This is general 
method which allows to generate ensemble of n-body configu- 
rations according to some prescribed probability density. That 
is, the probability of a given configuration in the ensemble is 
precisely that given by the probability density used to gener- 
ate "successive" configurations. In |0] this technique was used 
to modify directions of momentum vectors of number of se- 
lected particles from a system of n identical particles in order 
to impose the n-particle distributions derived from BE corre- 
lation functions. This procedure is then repeated many times, 
changing selected particles, until a kind of "equilibrium" is 
achieved. As shown in J2l one was able in this way to gener- 
ate typical multipion events which explicitly exhibit all corre- 
lations induced by Bose statistics. As a result of application 
of this algorithm a number of objects (called speckles) being 
clusters of large number of identical pions in the phase-space 
is being formed. The drawback of this method is that sym- 
metrization of clusters with sizes (represented by the num- 
ber of particles inside cluster, n c [ uster ) larger than n c i llster ps 10 
takes prohibitively long time. In |3J this method was there- 
fore modified by limiting symmetrization only to particles 
in clusters. This was possible by using wave packets to de- 
scribe produced particles instead of plane waves used in |Ql 
allowing therefore for localization of particles within certain 
phase-space volume and for providing the suitable criterium 
for defining a cluster. 

The accepting-rejection method investigated in [4] is based 
on the well known "hit-or-miss" technique of generating a 
set of random numbers according to a prescribed distribution. 
The method was designed to model collapse of a multiparticle 
wave function into a properly symmetrized state, as required 
by Bose quantum statistics. In contrast to the previous one it 
is sequential because n multiparticle event is constructed by 
first choosing single particle in phase space, then adding to 
it the second one according to the assumed 2-particle corre- 
lation function C2 x P2, then adding 3 th particle according to 
C3 P3 and so on. It is easy to realize that in this way one 
gets in the allowed phase space a "cell-like" structure because 
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regions with some particles inside them already present will 
have bigger chance to attract new particle. In a sense it looks 
like follows: first particle forms a seed for a first cell. When 
second particle is added to event it can, depending on its dis- 
tance from the first one, either remain in that cell or later on 
attract new particles and in this way start to form a new cell. 
This will then continue until all particles are used. Unfortu- 
nately, this sequential procedure is even more time consuming 
than the previous one. 



H. OUR PROPOSITION 

The above discussion shows that complexity of numerical 
symmetrization of wave function for all identical pions pro- 
duced in a given event (cf., eq.lfl])) can be substantially re- 
duced if only one can justify the idea that such symmetrization 
should be applied to groups of limited number of particles, 
as proposed in [3]. As argued there this can be achieved by 
dividing, according to some prescribed procedure, the initial 
set of secondaries into clusters consisting of identical parti- 
cles with similar momenta. Such groups of particles, but this 
time of (almost) equal energies, were introduced in [5] under 
the name Elementary Emitting Cells (EEC's). The physical 
justification was that bosonic particles, because of their sta- 
tistical properties, tend to occur, as much as possible, in the 
same state. One should realize now that such decomposition 
corresponds to factorization of permanent given by eq.© into 
matrix with a block structure: 



EECi 
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This idea of EEC was exploited by us further in [6] where 
instead of symmetrization of multiparticle wave function (de- 
pending on space-time positions, x, and energy-momenta, p) 
we worked in the number of particles basis. The bosonic char- 
acter of the produced secondaries means in this case the spe- 
cific bunching of identical particles in the phase space. In fact 
it is nothing else but modelling the correlations of fluctuations 
due to quantum statistics present in the system, which for 2- 
particle case are represented by: 
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o(n) is dispersion of the multiplicity distribution P(n) and p 
is the correlation coefficient depending on the type of parti- 
cles produced: p = + 1 , — 1 , for bosons, fermions and Boltz- 
mann statistics, respectively). The important feature of the 
EEC method is that number of particles in each cell follows 
by definition the geometrical (or Bose-Einstein) distribution 
characterized for identical bosonic particles. 

It should be mentioned at this point that the importance of 
bunching of particles in modelling quantum statistical effects 
has been demonstrated already in H. In this paper, follow- 
ing ideas of information theory, authors constructed Monte- 
Carlo (MC) event generator for multiparticle production pro- 
cesses (and applied it to e + e~ reactions). The main point 



was the assumption that particles of the same charge are lo- 
cated in cells (in their case they were constructed in rapidity 
space and were of equal size). It turned out that such divi- 
sion of phase space into cells was crucial for obtaining the 
characteristic form of the 2-body BEC correlation function 
Ci{Q~ \Pi~Pj\) (i.e., the one peaked and greater than unity 
at Q = and then decreasing in a characteristic way towards 
C% = 1 for large values of Q, out of which one usually deduces 
the spatio-temporal characteristics of the hadronization source 
Id). 

In our case we argue that the method of production of EEC 
proposed in [5] can be used to define the structure of clusters 
obtained in [3] (with, as it turns out, about n c i = 2 particles per 
cluster on average, depending on circumstances H). There- 
fore, instead of symmetrizing all particles in a given event, one 
symmetrizes separately particles in a number of EEC's with 
n c i <C n%< in tne plane wave approximation used in |01, one 
has therefore for some typical EEC the following « c /-particle 
probability function: 
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Notice that eq. (Q still contains (unmeasurable) positions of 
production of particles, {r a ( ; j}. Therefore they will be later 
eliminated by selecting them from some assumed distribution 
in a kind of numerical integration process. It corresponds to 
analytical integrations encountered before but in our approach 
we do not limit in any way the form of density function used 
(it had to be facto rizable before). It should be stressed that 
when all particles in the cluster are exactly in the same state 
then one gets, as expected, 
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However, in practice, even now from time to time one en- 
counters EEC with n c i ^s> 1, in which case eq.© is still very 
time consuming. For such cases we have to use some approx- 
imate schemes. The first one is sequential (cf. FIG. H}: one 
starts with some single particle in EEC and adds new particles 
one by one to all others already present and correlate them by 
using standard 2-particle probability: 1 + cos (Sr ■ hp). The 
resulting probability to get n c /-particle state is given by: 
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FIG. 1: Representation of algorithm for the sequential model, eq. 
©. 

Although appealing (in fact it resembles procedure used in 
yj]) it has drawback that when all particles are in the same 
state then 
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FIG. 2: Example of C2 modelled for single EEC with N = 4 by using 
2-particles (sequential) correlations as given by eq. $9$ (circles) and 
full A/-particles correlations as given by eq. 10 (stars). For compari- 
son, the analytical result, C 2 (SP) = 1 + sin 2 (R&P) / '{R&P) 2 , obtained 
for uniform one dimensional source in space with R = 1 fm is also 
shown (dotted-dashed line). 
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FIG. 3: Comparison of Ci modelled by using MC event generators 
with EECs (circles) and without EECs but with particles selected 
directly from the corresponding Bose-Einstein distribution |9|]. In 
both cases as reference event the Boltzmann distribution was used. 

i.e., for n c i > 2 the correlations are stronger and the maximal 
value of probability function gets bigger than allowed limit 
defined by eq.®. Both, exact (0 and sequential (O methods 
are compared in FIG. |2] Approximation (O keeps the width 
of C2 the same, however it differs substantially in C2(5P = 0) 
(from which one tends to estimate the so called chaoticity of 
the hadronizing source fljj] ) . The other possible approxima- 
tion, which was used in our hitherto applications [8], is that 



all particles in a given EEC are correlated only with the first 
particle defining this cell, not between themselves. It is inter- 
esting to notice that the results for C2 obtained this way are 
almost the same as those obtained by using the full method, 
eq. (0, but differ from that of eq. ((9). 

We close this section by noticing recent attempt to imitate 
the bosonic nature of particles produced by MC event gener- 
ator, in which they were choosen directly according to Bose- 
Einstein distribution [9]. This method seems to be very nat- 
ural and occurred to be also very fast. However, as one can 
see from FIG. [3] where it is compared with our approach, it 
leads to completely different C2 substantially underestimating 
the correlation function C2. We claim therefore that such pro- 
cedure inevitably loses some piece of important information, 
namely the fact that EEC's are formed and that Bose-Einstein 
distributions are only for particles in such EEC's, not in the 
whole event. In fact, signal of possible Bose-Einstein correla- 
tions without cells is seen in FIG.|3]are the trivial correlations, 
which can be eliminated by the proper choice of the reference 
event. 



m. SUMMARY 

We would like to summarize by stressing that there is no 
way to add to any of the existing MC event generators effects 
of quantum statistics, in particular Bose-Einstein one (BE). 
This is because they are all build on basis of classical physics 
with both the space-time and energy-momentum characteris- 
tics of produced secondaries used simultaneously. The only 
way out advocated here (albeit, most probably, not very prac- 
tical one and therefore hardly to be followed) is to build the 
multiparticle MC event generator ab initio, with BE properties 
(like bunching in phase space) being one of its basic principles 
and consisting its first step. All other features of such genera- 
tor would have to be added only after this. So far there is only 
one working example of this type [7], our efforts [8] aim for 
its further generalization and will be continued. 
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